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Abstract 

We investigate the pure penguin decays B — > 7T0 in the Standard Model(SM) 
and in the Constrained Minimal Supersymmetric Standard Model (CMSSM) using 
the QCD factorization approach and consider the Sudakov effects in the twist-3 
contribution. We find Br(B~ -> n'(j>) = (1.95 - 5.70) x 10~ 9 in SM and (1.1 - 



o 



2.4) x 10 in CMSSM with large tan/3 which is about one order of magnitude 



■ larger than that in SM. 
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1 Introduction 



One of charmless two-body nonleptonic decays of B mesons, the process B — > 07r, is 
interesting because it is a pure penguin process and, in particular, there are no anni- 
hilation diagram contributions the importance of which is still in dispute [M E]. It is 
sensitive to new physics due to all contributions arising from the penguin diagrams. The 
calculation of the hadronic matrix element relevant to the process is relatively reliable 
because of no contributions coming from diagrams of annihilation topology. Therefore, 
we shall investigate the process in both SM and MSSM. 

The study of exclusive processes with large momentum transfer in the perturbative 
QCD (PQCD) has been extensively carried out and it is shown that the application 
of PQCD to them is successful ||. The key point to apply PQCD is to prove that the 
factorization, the separation of the short- distance dynamics and long-distance dynamics, 
can be performed for those processes. Recently, two groups, Li et al. ||, [| and BBNS 0, 
have made significant progress in calculating hadronic matrix elements of local operators 
relevant to charmless two-body nonleptonic decays of B mesons in the PQCD framework. 
In the letter we shall use BBNS's method ( QCD factorization ) to calculate the hadronic 
matrix element of operators relevant to the decay B — > <J)tt. 

The decay B — > (pir has been studied by several people 0, [| [7|, ||, [|. The naive 
factorization or BSW model |T(J is used in calculating the hadronic matrix elements 
in Refs.0, [J. The modified perturbative QCD approach ||, instead of using BSW 
model, is used in Ref.Q. In Ref.||, the QCD improved factorization (simply, QCD 
factorization) is used and only the leading twist contribution is included. The nu- 
merical result of the branching ratio (Br) in SM given in Ref.[§] is about an order of 
magnitude larger than that in Refs.[]5], ||. That is, "non-factorizable" contributions (in- 
cluding vertex, penguin, and hard spectator scattering corrections), which are the 0(a s ) 
corrections to the leading order result, dominate over the leading order result. In Ref.[[| 
the same QCD factorization is used and the twist-3 contributions are included. However, 
the Br given in Ref.|| is (3-8) xlO -10 which is of the same order of or even smaller than 
that in Refs.|5|, EJ. Considering these disagreements, it is necessary to do a calculation of 
Br in SM using the QCD factorization. We carry out such a calculation in SM first and 
then in Constrained MSSM. The difference between our calculation and that in |9[] is how 
to calculate the twist-3 contributions. The authors in Ref.|| follow BBNS's approach, 
i.e., to introduce a phenomenological parameter instead of the integral containing end 
point singularity |IT| . Indeed, the solution of the end point singularity in investigating 
form factors of mesons is known for a long time [EJ . That is, to retain quarks' transverse 
momenta in both the hard scattering kernels and distribution amplitudes of mesons and 
to include the Sudakov suppression [H| make the integral convergent and computable. 
Thus, there is no any phenomenological parameter introduced. In the letter we shall use 
the method to calculate the twist-3 effects. Our numerical result of the Br(B ± — > 7T0) in 
SM using QCD factorization approach(QCDF) is Br(£- -> = (1.95-5.70) x 1(T 9 . 
We have also calculated the Br in Constrained MSSM in order to see supersymmetric 
(SUSY) effects on the decay. The numerical result in Constrained MSSM with large 
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tan/3 is Br(B~ — > 7r~0) = (1.1 — 2.4) x 10~ 8 depending on the choice of some relevant 
parameters. 

The AB = 1 effective weak Hamiltonian in SM is given by 

^eff = % E k{ciQI + c 2 qi+ E a<5i + c 77 g 77 + c 8g g 89 )+h.c, (i) 

V 2 p=n,c ^ i=3,...,10 ' 

where X p = V p t)V* d , Q\ 2 are the left-handed current-current operators arising from W- 
boson exchange, Qz,...fi and Qi^xo are QCD and electroweak penguin operators, and 
Q77 and Qs g are the electromagnetic and chromomagnetic dipole operators, respectively. 



Their explicit expressions can be found in, e.g., Ref. |Tl|. Follow BBNS approach J2|, the 
hadronic matrix elements of local operators Qi at the leading order of the heavy quark 
expansion can be written as 

(7r(p)0(ff)|g,|S(p)> = F B ^(q 2 ) C dvT\v)^{v) 

Jo 

+ f dtdudvT IJ (Z,u,v)$ B (Z)$ n (u)$f(v) (2) 
Jo 

where 4>m (M=0,7r, B) are light-cone distribution amplitudes of the meson M, Tf and 
T/ are hard scattering kernels. 

2 a s order corrections of hadronic matrix elements 

The a s order hard scattering kernels in Eq.(§) can be obtained by calculating the dia- 
grams in Fig.(l). Substituting the kernels into Eq.(§), we get 

ex C 

{■n(f)\Qi\B) asOTder = -^-—Fi^lQi-xlB)^ (3) 
where Fi=F for i=4,10, (-F-12) for i=6,8 and otherwise with 



/• = —ism-'' 1 v ' ' ' ' " 

In Eq. 



F = _121 n A_i 8 + / ' + / " (4 ) 



fi = Jdxg{x)%{x), (5) 
1 — 2x 

g(x) — 3 lnx — i3n (6) 

1 — x 

is the contribution from the diagrams (a)-(d), the vertex corrections, and fl presents 
the contribution from the hard spectator scattering diagrams (e) and (f ) which is of real 
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Figure 1: Order a s corrections to the hard-scattering kernels T 1 
and T n . 



non-factorization contribution. If we ignore the transverse momenta of partons, like that 
in Eq.(l), 



f^ 1 oc J d^dudv 



+ 



U V 



m B £ u 2 



(7) 



where the distribution amplitudes of n,<f),B mesons can be found in Refs. |17|, [18|, |19[ . 
When u — > 0, the twist-3 contribution, the next term of Eq.(^), will lead to divergence 
due to the end-point singularity. Therefore, it is necessary to consider the transverse 
momentum kx effect and include the Sudakov suppression factors to eliminate the end- 
point singularity ||13||. After including the transverse momenta of partons, Eq.(^) changes 
into 



n 



oc 



X 



J d^dudv GpkxG? 2 k 1 T<i 2 k2T 

-um 4 B (j) B (£)(f) n (u)(f)\\(v) 

[£um B + (k T - k 1T ) 2 ][-uvm B + (k T - k 1T + k 2 x) 
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+ 



\212 



(8) 



[£um% + (k T - k 1T )][— uvm B + (k T - k 1T + k 2 T 

The ^ resummation of large logarithmic corrections to B, (ft and tt meson distribution 
amplitudes leads to the presence of the exponentials S^S^ and S w respectively [F2U 

rt d/i 



s B (t) 

S v (t) 



exp 



exp 



Jl/b pi 



exp 



Jl/b 2 yU 

dp 



(9) 



-s(uP-, b 1 ) - s((l - u)P-, h) - 2 f ^ 7 («,(/i 2 )) 

Jl/ln fJ, 

with the quark anomalous dimension 7 = —a s /ir. The variables b, bi, and 62, corre- 
sponding to the parton transverse momenta h?, &it> and k 2 T respectively, represent the 
transverse extents of the B, it and (f) mesons, respectively. The expression for the ex- 
ponent s is referred to [23|]. The above Sudakov exponentials decrease so fast 
in the large h region that the B — > 7T0 hard amplitudes remain sufficiently perturbative 
in the end-point region. By a straightforward calculating, we obtain the hard spectator 
scattering contribution 

f n 4tt 2 f n f B 

Jd> 



N F^(ml)m 



B 



d^dudv J bdbb 2 db 2 



x 



um B V B {L b)V w (u, b)V\\(v, b 2 )K (-i^/uvm B b 2 ) 
6{b 2 - b)I (\f&m B b)K (^um B b 2 ) + 6(b - b 2 )I (^/^m B b 2 )K (^um B b) 
-2uvjj n m B V B (^ fo) - ^ 7 ^' — V\\(v, b 2 ) — K-^—iy/mirriBh) 
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2iJmim B 



x 



6{b 2 - b)I (^J^m B b)K (\ftum B b 2 ) + 6(b - b 2 )I (\/^m B b 2 )K (^J^um B b) 



(10) 



where f n (f B ) is the pion (B) meson decay constant, m B the B meson mass, F+^ 7r ( 
the B —>■ 7r form factor at the momentum transfer mi, £ the light-cone momentum frac- 
tion of the spectator in the B meson, Ii are modified Bessel functions of order % 
and Vb-i 'Pn, V<i> are B, ir, (ft corrected meson amplitudes with the exponentials S B , 
and S n respectively p0[ . As noted in Refs. [jP^| , although the twist-3 contribution is 
power suppressed it is numerically comparable with the twist-2 contribution due to the 
chirally-enhanced factor m£(/i) / 1 m B {lj)[fh u {ii) + fhd(fj)]. From Eq.([TD|), we obtain that 
the twist-3 contribution to fu is numerically about the fourth of the twist-2 contribution 
to fu. It is worth to note that the numerical result of the twist-2 contribution to fu 
obtained from Eq. fllOD is almost completely the same as that obtained without including 
the Sudakov factor, as expected. 
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3 The branching ratio in SM 



The effective Hamiltonian (HD results in the following matrix element for the decay 

G 



^<t>\H e n\B) = -^ £ \ P ^4>\T P \B) 

V * p=u,c 



(11) 



where 



= a 3 (ir(f)) (db) V -A ® (ss) V -a 

+ a 5 (n(p) {db)v^ A ® (ss) v+A 

+ a 7 (TT(p) (db)v^ A <%> \e s {ss) v+/ 

+ a 9 (7T0) {db) v ^ A ® |e s (ss) y _, 



(12) 



The symbol <S> indicates that the matrix elements of the operators in T p are to be eval- 
uated in the factorized form (ir<f)\ji ® j%\B) = (n\ji\B) (0|j 2 |O). The 0(a s ) corrections, 
including the nonfactorizable corrections corresponding to the diagrams (e) and (f) in 
Fig.l, of hadronic matrix elements are, by definition, included in the coefficients a^. 
Collecting the results in the above section, we have 



a 3 (vr0) 
a 5 (7T0) 
a?(vr0) 
ag(vr0) 



C3 + Jv Ci + ^W CiF ' 



Cs 
N 



a s Cf 
4>lv 



C 7 + ^ + ^C 8 (-F 



-12), 

12) + ^Pi m (C l + 3C 2 ) 



Jy 47r Jy 9tt 



2) , 



(13) 
(14) 
(15) 
(16) 



which has form of A + Ba s . In order to keep our calculation consistent, we use LO 
Wilson coefficients which contribute to B and NLO Wilson coefficients which contribute 
to A. In Eqs. (0-|D C F = (N 2 - 1)/(2JV) ( iV = 3 is the number of colors ), and P e p m 
arises from electroweak penguin contributions, Fig. 0(g), and is given by 



pp 
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4 / duu(l — u) In 



m„ 



q 2 u(l 



u 



(17) 



From Eq.flTl]), the decay amplitude for B — >• 7r </> is 



V2A(B° -j 

V * p=u,c 



7T 



% + a 5 - - (a? + a l) 



(18) 
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Table 1: Wilson coefficients Cj in the NDR scheme. Input parameters are 
Ag^ = 0.225 GeV, m t (m t ) = 167 GeV, m b (m b ) = 4.2 GeV, M w = 80.4 GeV, 
a = 1/129, and sin 2 ^ = 0.23. 



NLO 


Ci 


c 2 


c 3 




/■"""if 


/Hf 

c 6 


I <r\ 

H = m b /2 


1.137 


-0.295 


0.021 


-0.051 


0.010 


-0.065 


LI / / ify 


1 081 

-L . UO J. 


— n 1 QD 


n 014 

VJ.VJ It: 




n doq 


— D 049 


p, = 2rrib 


1.045 


-0.113 


0.009 


-0.025 


0.007 


-0.027 




C 7 /a 


C 8 /a 


C 9 /a 


C w /a 






\l = m b /2 


-0.024 


0.096 


-1.325 


0.331 






H = m b 


-0.011 


0.060 


-1.254 


0.223 






H = 2m b 


0.011 


0.039 


-1.195 


0.144 






LO 


Ci 


c 2 


c 3 


c 4 


c 5 


c e 


fi = m b /2 


1.185 


-0.387 


0.018 


-0.038 


0.010 


-0.053 


(jl = m b 


1.117 


-0.268 


0.012 


-0.027 


0.008 


-0.034 


H = 2m b 


1.074 


-0.181 


0.008 


-0.019 


0.006 


-0.022 




C 7 /a 


C 8 /a 


Cg/a 


Cio/a 






fi = m b /2 


-0.012 


0.045 


-1.358 


0.418 


-0.364 


-0.169 


fi = m b 


-0.001 


0.029 


-1.276 


0.288 


-0.318 


-0.151 


\i = 2m b 


0.018 


0.019 


-1.212 


0.193 


-0.281 


-0.136 



The relevant Wilson coefficients in NDR scheme are showed in Table [3] |§ . For the other 
parameters, we use 



f B = 0.190 GeV, U = 0.131 Gev, 
U = 0-237 Gev, fj = 0.215 Gev, F_ 



B — >TT j 



0.30 



and Wolfensein parameters fitted by Ciuchini as |14 
A = 0.819, 



A = 0.224, 

p = p(l- A 2 /2) = 0.224 , fj = 77(1 - A 2 /2) = 0.317 . 



(19) 



(20) 



The numerical results of the branch ratio at different scales are showed in Table [| 



4 Br in Constrained MSSM 

It has been shown that the neutral Higgs bosons (NHBs) do make significant contri- 
butions to leptonic and semileptonic rare B decays in Constrained MSSM with large 
tan p [^, |26], |27j . For b — > dss, it is expected that the similar enhancement of Br will 
happen since the mass of the strange quark is the same order as or a little of larger than 
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Table 2: Branch ratios of B 



it cf) at scale m&/2, m b and 2 raj,. 





\i = m b /2 


H = 


m b 


H = 


2ra& 




NF 


QCDF 


NF 


QCDF 


NF 


QCDF 


Br in SM 


2.47 x 1(T 9 


1.95 x 1(T 9 


6.53 x l(r lu 


4.51 x 1(T 9 


2.52 x 1(T 9 


5.70 x 1(T 9 


Br in CMSSM 




2.4 x 10" 8 




1.6 x 10~ 8 




1.1 x 10- 8 



that of muon. In the section we will calculate the Br of B — > 7r0 in the large tan (5 case 
of the Constrained MSSM. 

In addition to Eq. ([]]), we have |24], |25[ 



njnew 

lie 



G p 



L eff 



-At) £ QQ. + h.c. 



v2 i=n,...,i6 

where Qu to Qi6, the neutral Higgs penguins operators, are given by 

Qu = (dfys+pj^vi^s-p > Q12 = (dibj)s+p J2i (,QjQi)s-p , 

Q13 = (db) s +p J2q (QQ)s+p , Qu = (dibj)s+p J2i (QjQi)s+p , 

Qis = da^(l + l5 )bJ2 g q^^ + l5)q, 

QlQ = (^(7^(1 +J5)bj ^qqjVfivi 1 +75)?i- 

Here (9i92)s±p — 9i(l ±75)92- Then Eq.(p7|) is extended to 

(t^W* |B) = ^ E A P MT P + T p ncu |i?} , 

V ^ p=u,c 

where the term T™ eu arises from the neutral Higgs contributions, given by 
T p neu = r*(ji) [a n (7T0) + a Vi (Tx4>)}(db) v+A ® (ss) V - A ■ 



(21) 



(22) 



(23) 



(24) 



The Wilson coefficients Cq. (i=ll, 16) in Eq.(^TJ) is calculated in the same way as 
that in Ref.pBl, 27 j and results are 



C Qll (Mw) 



Cq 13 (M w ) 



CoAmw) 



Q,, 



m b m s tan 3 /3 2 ° . 

2^ U i2 1 UL K mb 



Ait Asm 2 e w M^X ti=lk=1 



■V2V i i(T VL K)* ka + Vt. 



(T UR fh u K) 



ks 



M w sin (3 



x (r hH + r A )Jx~^f B o (x x -,Xu k J + 0(tan 2 0), 

mfem^tan 3 /? ^ 6 

-a cm over4n 2 £ £ U i2 T UL K mb 

4 sin yu,i%A ( i=1 fc=1 

x (rhH - r A )^Jx~f B o (x^^Xu^ + 0(tan 2 /3), 

0, z = 12,14, 15, 16, 



V2y il (T [ ; i x)L + ^ 



(25) 

M w sin /? 

(26) 
(27) 
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where the definitions of various symbols are the same as those in Ref . . In calculating 
Cq^mw) the following values of relevant parameters are usedS 



tan/3 = 60, m h o = 110 Gev, m H o = 150 Gev, 

m H - = 200Gev, m b = 4.2 Gev, M 2 = 320 Gev, fi = 270 Gev, 

m tl = 120 Gev, 6~ ti = -tt/4, m s = 0.11 Gev. 

At the low scale we get 

C Ql2 + C Qli = 0. 139, 0.0897, 0.0565 for /i = m b /2, m b , 2m b . 

at in Eq. (||) in MSSM read as 

a n,i3 = —~r- -77 (fs + fs 1 ) Cqw.w > 



4vr N 



(28) 



(29) 



(30) 
(31) 



where 



// = -2 f 1 du[\n 2 u + 2\n u - 2Li 2 (u)} 4> s (u) 
Jo 



+2 du dk I db In 

JO JO 



\ 



4k, 



T + u 2 + u 



mi 



J (bk)V s (u,b) 



(32) 



j j _ 2-KVTLB f n f B 



[du}[db}S 2 (b 1 + b 2 )bV B (^b)V s (v,b 2 ) 



x [fi p (u + v) V p (u, h) + m B (£- v) V(u, 6i)] 



0(& 2 - b)I (bJu^m B ) K (bJu£,m B ) + 8(b - b 2 ) I (bJu£,m B ) K (\/-uvm b b 



(33) 



We found in numerical results that the last terms of the distribution amplitude <f> a (u) of 
meson make main contributions to an. 

The amplitude for for B~ — > 7r~</> now is given as 



A(B~ -> 7T" 
V * p=u,c 



«3 + «5 - ^ ( a 7 + a g) + r x ( a n + fl 13y 



where 



m B fi 



4e -p B 



x/^ J Pf- +w (mi)2^-p^4) 



(35) 



1 In CMSSM the mass spectrum and mixing of sparticles and Higgs bosons can be calculated given 
a set of values of a few parameters at the high (GUT or Planck) scale. Here, in stead of scanning the 
parameter space, we take reasonable values of relevant parameters for the sake of simplicity. 
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In calculating Oj (i=3,5,7,9) we should use the relevant Wilson coefficients in CMSSM. 
However, we still use their SM values in numerical calculations for simplicity because the 
contributions of SUSY only modify them in a few percents in most part of the parameter 
space including the values given above, Eq. (p8|) ||28|| . The numerical result of the Br is 
shown in Table ||. 



5 Summary 

In summary, we have studied the pure penguin process B~ — > n~(f) using QCD factoriza- 
tion approach, in particular, calculated the twist-3 contribution by including the Sudakov 
effects. We find Br(B~ — > 7r~</>) = (1.95 — 5.70) x in SM, which is roughly in agree- 
ment with that in Ref. |7[]. Comparing with the naive factorization (NF) result which is 
(0.7 — 3) x 10~ 9 , the QCD factorization result ( to the 0(a s ) ) is less sensitive to the decay 



scale, as can be seen from Table Bp. Actually, as noticed in Ref. [29], the coefficients a 



given in Eq. (|Tj-[T6|) are scale independent to the 0(a s ), which can be demonstrated by 
using the LO anomalous dimension matrix of relevant operators. However, from Table |2] 
one can see that there still is the significant scale dependence. Because we use the NLO 
Wilson coefficients the significant scale dependence comes mainly from the 0(a s ) cor- 
rections of hadronic matrix elements. Indeed the 0(a s ) corrections of hadronic matrix 
elements depend heavily on the scale (see, Eq.(§), which contains the factor 12 lnf/z/mft]). 
In order to decrease the scale dependence it is expected to calculate the a 2 s order cor- 
rections of hadronic matrix elements. We have also calculated the Br in Constrained 
MSSM with large tan/5 and the result is Br(£- -> 7r~0) = (1.1 - 2.4) x 10~ 8 . That is, 
the Br can be enhanced by one order of magnitude at most compared with that in SM, 
which is still far below the Babar experimental bound B r (B 0,± — > n 0,± <p) < 5.6 x 10~ 7 
at 90% CL M . 
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